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PROPERTIES OF CRYSTALLINE BISMUTH SELENIDE AND 
ITS USE AS A HALL EFFECT MAGNETOMETER 
by John A. Woollam, Harry Beale* and Ian L. Spain t 
Lewis Research Center 

SUMMARY 

Single crystals of both n- and p-type BigSeg grown by the Bridgman technique are 

found to make excellent Hall effect magnetometers. Plots of Hall resistivity p against 

yx 

magnetic field B to 10 tesla deviate from linearity by less than 1 percent and show no 
Shunbnikov-de Haas oscillations below 10 tesla at any temperature. The slope of the p 

yx 

against B curve varies by about 1 percent in the region 1. 1 to 35 K and by less than 20 
percent in the region 1. 1 to 300 K. Quenching tests demonstrate magnetometer reli- 
ability for cryogenic applications. 

From galvanomagnetic data the mobility and the carrier density as a function of tem- 
perature are found. Optical reflectivity experiments suggest a band gap of ~0. 08 eV. 

Analysis of all data indicates that the samples are nearly compensated (requiring at least 

25 

two -band conduction) and have semimetallic carrier density of ~10 per cubic meter. 

O 

Mobilities on the order of 10 square meter per volt per second at 300 K increase by 
greater than a factor of three at 1. 1 K. 


INTRODUCTION 

NASA currently has operating and under construction several superconducting mag- 
nets (ref. 1) to be used for plasma physics studies and other research related to space 
power and propulsion applications. Accurate and convenient measurement of high mag- 
netic fields becomes more difficult as higher field magnets and superconducting solenoids 
are developed. Nuclear magnetic resonance is very accurate, but high homogeneity of 
the magnetic field is needed (ref. 2), and it is often difficult to include a nuclear reso- 
nance measurement at the same time that other measurements are performed (ref. 3). 
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Copper magnetoresistors have nearly linear output with field at high fields, but their cal- 
ibration is very temperature sensitive and a quadratic component is present in the output 
at low field (ref. 4). Hall effect devices are small (dimensions on the order of a few mm) 
and quite linear in field at high temperature. The best Hall devices available up to now 
are films of indium arsenide (ref. 5), but these are frequently unreliable upon cycling to 
low temperature and often exhibit Shubnikov-de Haas oscillations at low temperatures 
(ref. 6). 

We have found single crystals of BigSe^ grown by the Bridgman technique to be ideal 
magnetometers in magnetic fields to at least 10 tesla and over a wide range of temper- 
ature. hi this report, the measured Hall resistivity is reported for a representative p- 
type sample as a function of temperature from 1. 2 to 300 K and as a function of magnetic 
field to 10 tesla. Magneto resistivity was also measured and is used with the Hall effect 
results to help characterize the magnetometer material. A series of direct immersions 
from 300 K into cryogenic liquids demonstrated reliability for low temperature uses. 

In addition to its use as a magnetometer, which we have demonstrated in this report, 
B^Seg is of practical interest for several other reasons. The thermoelectric figure of 

e^ making these materials useful for ther- 
moelectric energy conversion. Secondly, the band gap is found near 0. 08 eV in our sam- 
ples. Thus, this material is of possible interest for infrared radiation detectors (ref. 8). 
The optimization of the material for these applications will require further study and 
improvement of the electronic properties. This report describes studies of electronic 
and optical properties of n- and p-type single crystals in high magnetic fields and over a 
wide range of temperatures. Hashimoto (ref. 9), and Caywood and Miller (ref. TO) have 
studied galvanomagnetic properties, but their studies were limited to temperatures above 
80 K and to low magnetic fields. Optical studies were made by Greenaway and Harbeke 
(ref. 11) who find an energy band gap of =0. 16 eV in their samples. 


merit is high (ref. 7) for both B^Se-j and B^T 


MATERIALS PREPARATION AND MEASUREMENT TECHNIQUE 

The starting materials for this series of experiments were 99. 999 percent pure bis- 
muth and selenium shot. To make p-type samples, 100 grams of material were weighed 
out, to an accuracy of ±0. 0005 gram, for a composition of 60. 25 atomic percent selenium 

and 39. 75 atomic percent bismuth. The mixture of shot was placed in a vycor ampoule, 
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evacuated to approximately 10 millimeters of mercury and sealed. The ampoule was 
maintained at 850° C (the melting temperature is 710° C) for 12 hours in a furnace vibrat- 
ing at approximately 8 hertz. Samples were made using a Bridgman furnace (ref. 12), 
with a growth rate of 0. 012 meter per hour. The product was a polycrystallirte boule, 

0. 07 meter long by 0. 015 meter diameter. Grains were large in cross section and typ- 
ically the length of the boule. By chemical analysis, the composition varied from 39. 97 
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atomic percent bismuth near the bottom of the boule to 39. 82 atomic percent near the top. 
The boule was sectioned into 0. 018 meter lengths using a wire saw with a diamond im- 
pregnated wire. One of these shorter cylinders from the upper end of the boule was 
etched, using a mixture of methanol and bromine, to reveal the grain boundaries and then 
cleaved to produce single crystals. The predominant cleavage plane was the plane per- 
pendicular to the trigonal axis. 

Specimens were cut by sand erosion to the shape of a conventional Hall resistivity 
sample as shown in figure 1. Typical sample dimensions were 0. 05 by 0. 14 by 
0. 80X10 meter. BigSe^ has a trigonal unit cell and exhibits the R3m symmetry. If the 
z-axis is taken as the trigonal axis and the x-axis as the long axis of the sample, then 
the specimens were oriented and cut so the xz -plane was one of the three equivalent re- 
flection planes. The samples were again etched and re-examined for cracks between the 
layer planes using a simple light microscope. Leads were attached with a conductive 
epoxy, carefully covering the entire tip region. This procedure eliminated any spurious 
effects due to the pronounced layer structure of the crystal. The previously described 
procedure produced p-type material as grown; n-type was produced in a second boule 
nominally 0. 05 percent bismuth rich by annealing for 48 hours at 600° C. 



Figure 1. - Sample geometry showing tabs for attachment of electrical leads 
with conductive epoxy. (Electroplated copper was also successful. ) 
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A dc constant current supply, having an output calibrated to one part in 10 , pro- 
vided currents of typically 100 milliamperes. Voltages were measured by using a dc 

5 

amplifier with gains accurate to one part in 10 and were plotted on a calibrated x-y re- 
corder. A water-cooled solenoid, calibrated against a rotating coil magnetometer and a 
standard magnet, provided the magnetic field. Field values ranged between zero and ±10 
tesla (1 tesla = 10 kG), with the field oriented along the trigonal axis of the sample. 
Temperatures between 1. 1 and 300 K were achieved by using liquid helium and liquid ni- 
trogen with a heater -controlled system employing a gallium arsenide temperature sensor 
(ref. 13). 
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THEORETICAL CONSIDERATIONS 


The experimentally measured quantities were the Hall resistivity p and the mag- 

xy — 

netoresistivity p . For the magnetic field along the trigonal axis of the R3m structure 

AA 

of Bi„Se„, the conductivity tensor elements or and or can be obtained from the fol- 
di o xx xy 

lowing relations: 
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(Symbols are defined in appendix A.) These conductivity tensor elements can be ex- 
panded (ref. 14) as 
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(3) 




and 
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(4) 


where the sum is over the different carriers and where n. is the concentration of the 

1U 11* 1 

i ul carrier, e^ is the charge of the i Lr carrier, B is the magnetic flux density, and 
Qf| and are degeneracy parameters, the values of which do not usually deviate from 
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unity by greater than 20 percent. For the remaining development of this section, these 
parameters will be taken as unity, as a simplification. By definition, 



(5) 


where ju^ is the mobility of the i** 1 carrier. 

The simplest analysis allows the use of a one carrier model. In this case, only two 
quantities are necessary to specify the system, the carrier concentration and charge 
mobility. If 


Ap 3 Pxx^ ~ ^xx^ ^ 

Po " P*x< B - °> 


then the magnetoconductivity mobility is 


( 6 ) 


Pc 



m^/(V)(sec) 


(7) 


where B is in tesla and V is in mks practical volts. Equation (7) follows from equa- 
tion (4) for a one -band system assuming p » p . The latter was experimentally 

xx xy 

found to be true in BigSe^. However, as will be shown, BigSe^ can be described as a 
two -band system. The carrier concentration is 


n = 


<W B - 



( 8 ) 


which follows from equation (4) if B « Bj, which was also found experimentally in 
BigSeg. The system is then determined. However, another independent operational rela- 
tion is available from the Hall parameter. With 

Rjj(B) = ^ (9) 


the Hall mobility is 


p R = " 0) 


( 10 ) 
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Therefore three relations (eqs. (7), (8), and (10)) exist for the determination of two quan- 
tities (n and ju). These relations can be used to assess the applicability of the chosen 
model (either single carrier or multicarrier). If the quantities /i c and jUjj differ by 
less than say 20 percent, then it may be inferred that the single carrier model is appli- 
cable to the experimental data. Such a small discrepancy implies the degeneracy param- 
eters a and /3 deviate from unity. A small discrepancy is also expected dependent on 
the predominant scattering mechanism. However, if the disparity between the fi and 
jijj is larger, then it may be inferred that more than one carrier is responsible for the 
measured galvanomagnetic properties. This results from equations (7) and (10) being 
valid for the one carrier case only. 

Consider a two carrier model. If the degeneracy parameters are again taken as 
unity, then only four quantities are necessary to specify the system, instead of eight. 

With the definitions 



and 
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H = ^2 

H 2 B x 
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( 12 ) 


the two carrier analogs of equations (7), (8), and (9) are, from Soule (ref. 15), 


and 


M, 


VWW (p o B 2 < 
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x 10 4 m 2 /(V)(sec) 


a xx (B = 0)=n ll e l^l( 1 + g) 


(13) 


(14) 


1^(0) = J- < a - b J. 


n l e (a + b) 2 


(15) 


for B « B^. For equation (15), e 2 = _e i • Therefore, only three relations exist (eqs. 
(13), (14), and (15)) for the determination of four quantities (jz-, \i 2 , n^, and ng). An 
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additional relation can be obtained from the high field dependence of the Hall parameter. 
This may be accomplished by fitting to the magnetic field dependence of Rjj(B) if data is 
only available in the range B < Bjj or if data to the high field limit (B » B^) is avail- 
able, then with the aid of the formula 


R (°°) = r-r-^ 

|e|(n 2 - nj) 


(16) 


Equations (13) to (16) can then be used to determine the four unknowns. 


RESULTS 

The data discussed in this report are on p-type material. We also have studied n- 
type materials, and the results are very similar. All of the discussion concerning mag- 
netometer applications is equally valid for both n- and p-type. 

Representative plots of the Hall resistivity as a function of magnetic field B for 
various temperatures are shown in figure 2 for the p-type material. Over the range of 

field from zero to 8. 0 tesla, the variation of p is linear in B. Figure 3 displays the 

yx 

variation in p-type material of p with temperature at a fixed field of 7. 0 tesla. This 

yx 

is essentially a plot of the temperature dependence of the magnetometer sensitivity. In 
table I are listed the values of the Hall and magnetoconductivity mobilities which were 
computed using equations (7) and (10). Clearly these quantities differ significantly, indi- 
cating the presence of more than one carrier. 

Figure 4 is a composite, displaying the temperature dependence for p-type material 
of p ( o) and also the magnetic field dependence of p (B) for various temperatures. 

The weak field dependence of p is very evident, being ~15 percent at 4.2 K and ~2 

AA 

percent at 300 K at 7 tesla. It is also evident that there is a quadratic component of field 

dependence. From figures 3 and 4, note that p is always less than 0. 8 percent of 

xy 

p . The resistivity ratio p (300 K)/p (4. 2 K) = 4 for this sample. Equation (6) was 

'XX XX XX 

used to compute the magnetoresistance which is shown in figure 5. The magnetoconduc- 
tivity mobility was computed using equation (7) and is also shown in figure 5. An average 
carrier concentration was computed using equation (8) and is shown in figure 6. 

Finally, we have measured the reflectivity as a function of wavelength at 300 K. A 
small rise in reflectivity near 15 micrometers suggests the presence of a band gap of 
0. 08 eV, which is considerably smaller than the value of 0. 16 eV found by Greanaway 
and Harbeke (ref. 11). 
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TABLE I 4 - HALL MOBILITIES, CONDUCTIVITY MOBILITIES, CON- 
DUCTIVITIES, AND TOTAL CARRIER DENSITY FOR A REPRE- 


SENTATIVE p-TYPE SAMPLE FROM AN UNANNEALED BOULE 


Temperature, 

K 

Hall 

mobility, 

, y H- 

m V(V)(sec) 

Magnetoconduc - 
tivity mobility, 

2 Mc ’ 

m V(V)(sec) 

Conductivity 
at zero field, 

CT xx (0) > 

(ftm) -1 

Net carrier 
density, 

n l + *12’ 

4.2 

1.3X10" 2 

4. 75X10" 2 

1. 11X10 5 

1.5X10 25 

38 

1. 3 

4.40 

1. 11 

1.6 

78 

.88 

3.50 

.714 

1.3 

105 

.79 

3. 10 

.589 

1.2 

155 

.60 

2.40 

.455 

1.2 

198 

— 

2.3 

.385 

1 . 1 

238 

— 

2. 1 

.345 

1.03 

300 

.42 

1.75 

.286 

1.00 
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Figure 4. - Resistivity plotted against magnetic field for various temperatures and 
against temperature. 






DISCUSSION OF RESULTS 


Applications as a Magnetometer 

The most significant results of the present investigation for practical applications, 

is that the Hall resistivity p is linear in B, within 1 percent, over the temperature 

yx 

range 1. 1 to 300 K. The data for several temperatures are conveniently displayed in 
figure 3. 

The linear Hall resistivity with magnetic field may be explained using the theory dis- 
cussed earlier, as will now be shown. The result of inverting equations (1) and (2) is 



and 



These equations may be expanded by using equations (3) and (4) as (see appendix B for 
details) : 



and 
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where 


M i sl+ (~j (21) 

Consequently, whether p x y has strong linear or cubic field dependence is controlled by 
the quantities (1 - a) and (b 2 - a). As "a" approaches unity (meaning the material is al- 
most completely compensated, n^ « ng), p xy will become nearly linear in field. The 
parabolic behavior of p is demonstrated by equation (20). Again, linearity of p 
with B and a quadratic component of p with B were observed experimentally. The 
following experimental results will be used to evaluate Bj, Bg, and a in equations (19) 
and (20): 


Rjj(B) - Rjj(0) 
RjjfO) 


B=7 tesla 
T=4.2 K 


~ 0.015 


Pxx (B) - p„(0) 


XX' 


Pxx (0) 


B=7 tesla 
T=4. 2 K 


~0. 15 


( 22 ) 


(23) 


Pxx (B) 

— — > 0.008 (24) 

Pxy( B > 

From the measured values of mobility (table I) the B. ’s are large compared with the 
largest B used (see eq. (5)). Thus equations (19) and (20) yield 


Rgg* - Kg* 0 * B 2 (l-a) 


Rh(0) 


B 2 (b 2 - a) 


Pxx^ ~ p xx^ 1 / B \ 2 (ab+1) 


Pxx (0) 


b\B ly / (a + b) 


Pxx (B) 

*>xy( B > 


b(a + b) + flL) (ab + 1) 


1 / 

B(b 2 - a) 


B„ 


(25) 


(26) 


( 27 ) 
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There are then three equations in three unknowns, a, Bj, and Bg which could be solved. 
To simplify solution, various trial values of a and b were chosen and varied until equa- 
tions (25), (26), and (27) fit results (22), (23), and (24) to within a few percent. Incom- 
patibility of solutions for a « 1 was found. But with a s 1, the following results were 
found for a temperature of 4. 2 K: 


a 0.99974 
Bj 18. 00 tesla 

> 

Bg 18.02 tesla 
b 1.001184 


(28) 


Thus the linearity of p with B (see eq. (19)) results principally from the near com- 

xy 

pensation of carriers (n^ “ ng or a s 1). By using equations (28) and (5), the mobilities 
should then be 


jlij_ “ p 2 “ m 2 /(V)(sec) 


(29) 


Experimental values of p for various temperatures are listed in table I. 

Nearly as significant as the linearity of p with B and the absence of the 

yx 

Shubnikov-de Haas effect below 10 tesla in BigSe^ is the very small temperature variation 
of the slope of p , _ against B. This means that a magnetometer needs only to be cal- 
ibrated at a few temperatures (e.g. , 4.2, 78, and 300 K) and can be used for accurate 
field measurements over wide ranges of temperature. This property contrasts with the 
strong temperature dependence of some presently used Hall effect magnetometers 
(ref. 5), and all copper magnetoresistance magnetometers (ref. 4). Copper magneto - 
resistance typically changes by 100 percent over the range of 35 to 4. 2 K. 

A Hall effect magnetometer has an inherent advantage over magnetoresistance mag- 

o 

netometers in that the magnetoresistance of most materials is proportional to B at low 
B. This is true, for example, in copper magnetometers which are fairly linear only 
above a magnetic field of typically 0. 2 to 1.0 tesla (ref. 4). 

The most serious problem with previous Hall magnetometers was the change in 
characteristics, or complete failure, after cycling to cryogenic temperatures. For 
this reason BigSe^ probes have been cycled in the following manner: 

(1) Ten direct immersions into liquid nitrogen from 300 K with the Hall resistivity 
measured to ±5 tesla each time the probe was in liquid. All values of p^ at 5 tesla 
were within 1 percent of each other. 

(2) Ten direct immersions into liquid helium from 300 K with the Hall resistivity 
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measured to ±5 tesla each time the probe was in liquid. All values of p^ at 5 testla 
were within 0. 2 percent of each other. 

(3) Thirty direct immersions into liquid nitrogen from 300 K with p^ measured 
in liquid helium at 5 tesla after each 10 immersions into liquid nitrogen. All values of 
p at 5 testla were within 0. 3 percent of each other. 


Characterization of the Material 

In the last section it was demonstrated that the data can be interpreted by a two -band 
model where a = 0.999, = Bg = 18 tesla, and mobilities are ~0.056 square meter 

per volt per second. Figure 4 shows that the total resistivity decreases by a factor of 
four on going from 300 to 1. 2 K, and figure 5 shows that mobility increases by a factor of 
about three as temperature decreases. Elementary two -band conduction yields 

ff xx =n l e * i l +n 2 e * i 2 (30) 

Since = cr^ in our samples, there is a slight increase in the number of carriers as 
temperature decreases as is shown in figure 6. This behavior is expected for a degen- 
erate material. 

Since the optical reflectivity suggests a band gap of ~0. 08 eV, a band model shown in 
figure 7 is proposed. A degenerate n-type impurity band overlaps the conduction band. 



Figure 7. - Density of states model used to qualita- 
tively explain data. 
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These bands are located at different positions in k space, and carrier compensation (as 
in a semimetal) can thus exist. From our analysis, the mobilities of the two bands are 
nearly the same. A band gap exists between the impurity band and the conduction band. 
Changes in optical reflectivity are then caused by optical excitation of carriers to the con- 
duction band. Transmission experiments were not possible due to strong absorption. 

The number of holes and the number of electrons could not accurately be found from 
the two-band model since high enough magnetic fields (> 18 tesla) were not available to 
make equation (16) applicable. An estimate of n^ + ng was made, however, using equa- 
tion (30) from the two -band model, and these values are listed in table I as a function of 

2*5 -3 

temperature. Then, knowing a s 0.9997, n^ s ng = 10 meter . These numbers are 
characteristic of semimetallic densities and are consistent with the model of a heavily 
doped semiconductor proposed. Since chemical analysis revealed very low concentrations 
of anything but bismuth and selenium, the impurity states are created by lack of 
stoichiometry. 


CONCLUSIONS 

It has been found that BigSe^ single crystals grown by the Bridgman technique are ex- 
cellent Hall effect magnetometers for the temperature range 1. 1 to 300 K and magnetic 

fields to at least 10 tesla. The slope of the .p against B curve varies by less than 

yx 

1. 0 percent in the temperature range 4 to 35 K. In the same temperature range, a copper 
magnetometer calibration changes by 100 percent. Quenching tests demonstrate mag- 
netometer reliability for cryogenics applications. 

Our samples of BigSe^ are highly degenerate with an impurity band partially filling 
the band gap of the intrinsic material. This leaves a band gap of approximately 0.08 eV. 
Conduction is found to be by at least two bands, with mobilities of about 0.056 square 
meter per volt per second, and near compensation of carriers. 

To be useful as optical devices (lasers or radiation detectors, for example), crystals 
with greater stoichiometry need to be grown. This will reduce the free carrier optical 
absorption and increase the mobility. For possible laser use, the position of maxima and 
minima in the energy bands in k space must also be found. It is interesting that BigSe^ 
would not have been useful as a Hall effect magnetometer if the low mobility and lack of 
stoichiometry had not existed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 8, 1972, 

503-10. 
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APPENDIX A 


\ ■ 
v 


a 


B 

B i 

b 


c 

e 


n i 

n l +n 2 

Rh(B) 

M 

Mb 

M c 

M c 

Mi 

Mxx 

Pxx< B > 

Pyx 

Ap/p 

*xx< B > 


SYMBOLS 

ratio of concentrations of two carriers, ng/n^ 
magnetic field 

characteristic field of i** 1 carrier 

ratio of mobilities of two carriers, = B 2 / ^ B 1 

velocity of light 

electronic charge 

charge of i^* 1 carrier 

number of carriers per unit volume 

number of carriers per unit volume for i ul carrier 

net carrier density 

Hall constant, p / B 
xy 

mobility 
Hall mobility 

magnetoconductivity mobility 
average magnetoconductivity mobility 
mobility of i** 1 carrier 
resistivity 

resistivity as function of B 
Hall resistivity 
magnetoresistance 
conductivity at field B 
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APPENDIX B 


DERIVATION OF EQUATIONS (19) AND (20) 

Using equations (3) and (17) and definitions (11) and (12) gives 



where M. are defined by equation (21). Equation (32) then reduces to the desired 
expression 



It is easy (but lengthy) to show that if B « Bj or if a^b “1, then there will be no 
field dependence of the denominator. 
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Using equations (4) and (18) results in 
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